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We investigate small— a; QCD effects in forward-jet production in deep inelastic scattering in the 
kinematic regime where the virtuality of the photon and the transverse momentum of the jet are two 
hard scales of about the same magnitude. We show that the data from HERA published by the HI 
and ZEUS collaborations are well described by leading-logarithmic BFKL predictions. Parametriza- 
tions containing saturation effects expected to be relevant at higher energies also compare well to 
the present data. We extend our analysis to Mueller-Navelet jets at the LHC and discuss to what 
extent this observable could test these small— a; effects and help distinguishing between the different 
descriptions. 
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CN ■ I. INTRODUCTION 



The Regge limit of perturbative QCD comes about when the centre-of-mass energy in a collision is much bigger 
than the fixed hard scales of the problem. In this limit usually called the small— a; regime, parton densities inside the 
projectiles grow with increasing energy, leading to the growth of the scattering amplitudes. As long as the densities 
are not too high, this growth is described by the Balitsky-Fadin-Kuraev-Lipatov (BFKL) equation |l| that resums 
the leading logarithms. As the parton density becomes higher and the scattering amplitudes approach the unitarity 
limit, one enters a regime called saturation 0, El El 01 where the BFKL evolution breaks down and parton densities 
saturate. 
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In the past years, as colliders started to explore the small— a; regime, proposals were made to test the relevance of 



the BFKL equation at the available energies. In this paper we concentrate on two of the proposed measurements: 
forward jets [g in deep inelastic scattering (DIS) and Mueller-Navelet jets [tJ in hadron-hadron collisions. Forward-jet 
, production is a process in which the virtual photon interacts with the proton and a jet is detected in the forward 
direction of the proton. The virtuality of the photon and the squared transverse momentum of the jet are hard 
scales of about the same magnitude. In the case of Mueller-Navelet jets, a proton interacts with another proton or 
antiproton and a jet is detected in each of the two forward directions; the transverse momenta of the jets are as well 
hard scales of about the same magnitude. If the total energy in the photon-proton (for forward jets) or proton-proton 
(for Mueller-Navelet jets) collision is large enough, these processes feature the kinematics corresponding to the Regge 
limit. 

The description of forward jets with fixed-order perturbative QCD in the Bjorken limit amounts in the following. 
Large logarithms coming from the strong ordering between the soft proton scale and the hard forward-jet scale are 
resummed using the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution equation Q and the hard cross- 
section is computed at fixed order in the coupling constant. In the small— a; regime, due to the extra ordering between 
the total energy and the hard scales, other large logarithms arise and should be resummed within the hard cross-section 
itself. In other words, the inclusion of small— a; effects aims at improving QCD predictions by replacing fixed-order 
hard cross-sections with resummed hard cross-section, using the BFKL equation or, at even higher energies, using 
resummations that include saturation effects. 

To study different observables in this small— a; regime, a convenient approach is to formulate the cross-sections in 
terms of scattering amplitudes for colorless combinations of partons. The simplest of those is the qq dipole, a quark- 
antiquark pair in the color singlet state; it describes for instance the interaction of a virtual photon. Any colorless gg, 
qqg, ... multiplets can a priori be involved, for instance the gluon-gluon (gg) dipole is what describes gluon emissions. 

To compute the evolution of those scattering amplitudes with energy, the QCD dipole model Q has been developed. 
This formalism constructs the light-cone wavefunction of a qq dipole in the leading logarithmic approximation. As 
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the energy increases, the original dipole evolves and the wavefunction of this evolved dipole is described as a system 
of elementary qq dipoles. When this system of dipoles scatters on a target, the scattering amplitude has been shown 
to obey the BFKL equation. Interestingly enough, the dipole formalism was shown to be also well-suited to include 
density effects and non-linearities that lead to saturation and unitarization of the scattering amplitudes 0, ^| . This 
is why the dipole picture is suitable for investigating the small— £ regime of QCD, it allows to study both BFKL and 
saturation effects within the same theoretical framework. 

The formulation of the forward-jet and Mueller-Navelet jet processes in terms of dipole amplitudes has been ad- 
dressed in [l^ ■ I n both cases the problem is similar to the one of onium-onium scattering: the growth with energy of 
the total cross-section due to BFKL evolution is damped by saturation effects which arise purely perturbatively. For 
instance in the large— N c limit, that involves multiple Pomeron exchanges. In our study, we consider both the BFKL 
energy regime and saturation effects. We shall implement saturation in a very simple way, using a phenomenological 
parametrization inspired by the Golec-Biernat and Wiisthoff (GBW) approach [Lj which gave a good description of 
the proton structure functions with a very few parameters. 

First comparisons of small— a; predictions with forward-jet data from HERA were quite successful: the first sets of 
data published by the HI Q and ZEUS 0] collaborations are well described b y th e leading-logarithmic (LL) BFKL 
predictions ^(| and also show compatibility with saturation parametrizations [VA Il8| . In both cases, the DGLAP 
resummation associated accounted for leading logarithms. 

Very recently, new forward-jet experimental results have been published l2(lj|. They involve a broader range of 
observables with several differential cross-sections and go to smaller values of x than the previous measurements. As 
QCD at next-to-leading order (NLO) is not sufficient to describe the small— x data, we shall address the following 
issues: whether the BFKL-LL predictions keep being in good agreement and whether saturation parametrizations 
still show compatibility. The first part of the present work is devoted to those questions. 

In the second part of the paper, we deal with Mueller-Navelet jets. We display BFKL-LL predictions in the LHC 
energy range for different differential cross-sections. We compare them with saturation predictions obtained from our 
parametrizations of saturation effects constrained by the forward-jet data. We propose different measurements and 
discuss their potential for identifying BFKL and saturation behaviors. 

The paper will be organized as follows. In section II, we compute the forward-jet cross-section in the high-energy 
regime and express it in terms of a dipole-dipole cross-section. We compare the BFKL-LL predictions and the 
saturation parametrization with the new HI and ZEUS data for several differential cross-sections. In section II, we 
compute the Mueller-Navelet jet cross-section and show the BFKL and saturation predictions for LHC energies and 
for several differential cross-sections. The final section V is devoted to conclusion and outlook. 



II. FORWARD-JET PRODUCTION 

Forward-jet production in a lepton-proton collision is represented in Fig.l with the different kinematic variables. 
We denote yfs the total energy of the lepton-proton collision and Q 2 the virtuality of the intermediate photon 
that undergoes the hadronic interaction. We shall use the usual kinematic variables of deep inelastic scattering: 
x = Q 2 1 (Q 2 + W 2 ) and y = Q 2 /(xs) where W is the center-of-mass energy of the photon-proton collision. In addition, 
kx 2> Aqcd is the jet transverse momentum and xj its longitudinal momentum fraction with respect to the proton. 
In the following, we compute the forward-jet cross-section in the high-energy limit, recall the BFKL predictions and 
give our formulation of the saturation model. 

A. Formulation 

The QCD cross-section for forward-jet production reads 

d^v = aern J i f ' da^ JX | d*£^ \ | da^ JX y 2 \ 

dxdQ 2 dx j dk 2 -, nxQ 2 1 I dxjdk 2 -, dxjdk 2 ^ J dxjdk 2 -, 2 J ' 

where da^*^~* JX jdxjdk\ is the cross-section for forward-jet production in the collision of the transversely (T) or 
longitudinally (L) polarized virtual photon with the target proton. 

We now consider the high-energy regime i<1. In an appropriate frame called the dipole frame, the virtual photon 
undergoes the hadronic interaction via a fluctuation into a dipole. The dipole then interacts with the target proton 
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FIG. 1: Production of a forward jet in a proton-lepton collision. The kinematic variables of the problem are displayed. Q 2 is 
the virtuality of the photon that undergoes the hadronic interaction, s and W 2 are the total energies squared in the lepton- 
proton and photon-proton collisions respectively, kr is the transverse momentum of the forward jet and xj is its longitudinal 
momentum fraction with respect to the incident proton. Y is the rapidity interval between the two hard probes. 



and one has the following factorization 

The wavefunctions ipj* an d V'Z describe the splitting of the photon on the dipole and da q qjdxjdk\ is the cross-section 
for forward-jet production in the dipole-proton collision. i\hL and ipl are given by 



\r T (r, z; Q)| 2 = Yl 4^ + (1 " z?)z{l-z)Q 2 Kl(VW^)Q\r\) (3) 

/ 



WliT, z; Q)| 2 = 2=^£ £ ~ zfKl (v^(l=5)g|r|) (4) 

/ 

for a transversely © and longitudinally polarized photon where e/ is the charge of the quark 1 with flavor /. 
The integration variable r is the transverse size of the qq pair and z is the longitudinal momentum fraction of the 
antiquark with respect to the photon. In the leading logarithmic approximation we are interested in, the cross-section 
dtjqq/dxjdk^ does not depend on z but only on the dipole size r. This cross-section has been computed in 1 1 "2f where 
it was shown that the emission of the forward jet can be described through the interaction of an effective gluonic (gg) 
dipole: 

with Y = log(xj/x) the rapidity assumed to be very large. &(qq)(g g ) ( r > r 5 Y) is the cross-section in the collision of a 
qq dipole of size r with a gg dipole of size f with total rapidity Y. f e ff(xj,k^) is the effective parton distribution 



1 We consider massless quarks and sum over four flavors in J^J and J1J. This is justified considering the rather high values of the photon 
virtuality (Q 2 >5 GeV 2 ) used for the measurement. 



function and resums the leading logarithms log(k T / Aq C D ) . It has the following expression 

feff(xj,k T ) = g(x Jl k T ) + ^(q(xj,k T ) + q(xj,k T )) , (6) 

where g (resp. q, q) is the gluon (resp. quark, antiquark) distribution function in the incident proton. 

Let us comment formula (0- Since the forward jet measurement involves perturbative values of kx and moderate 
values of xj, it is not surprising that formula J^J) features the collinear factorization of feff, note also that k T has 
been chosen as the factorization scale. The remaining hard interaction is between a gg dipole and the incident qq 
dipole of size r. The gg dipole emerges as the effective degree of freedom for the gluon emission at high energies 12]. 
This feature has been pointed out several times [2l| . 

Formulae Q-© express the forward-jet observable (JIJ in terms of the cross-section (7(qq)(gg) which contains the 
high-energy QCD dynamics: the problem is analogous to the one of onium-onium scattering. In the next subsection 
(II-B), we deal with the BFKL energy regime for which the interaction between the qq dipole and gg dipole is restricted 
to a Pomeron exchange. In that case of course, our formulation is equivalent to the kx— factorization approach. In 
subsection (II-C), we go beyond kx— factorization and investigate the saturation regime in which &( q q){ gg ) a priori 
contains any number of gluon exchanges. 



B. The BFKL energy regime 

The BFKL qq-gg dipole-dipole cross-section reads (see for instance |2^ ) 

bfkl , - V \ o 2 2 f d l {r/rf 1 ( a s N c \ 

a (^)( 9a )^^ Y ) = 2 ™s r J ^- 7 2 (1 _ 7 )2 eX P { — XW Y J ( ? ) 

with the complex integral running along the imaginary axis from 1/2 — ioo to 1/2 + ioo and with the BFKL kernel 
given by 

X ( 7 ) = 2^(1) - V>(1 - 7) - V>(7) (8) 

where ip{"j) is the logarithmic derivative of the Gamma function. It comes about when the interaction between the 
qq— dipole and the gg— dipole is restricted to a two-gluon exchange. Summing the leading-logarithmic contributions 
of ladders with any number of real gluon emissions, one obtains the BFKL Pomeron and the resulting growth of the 
cross-section with rapidity. 

Inserting (J7J) in J5J) and (J2J), one obtains 

da T *r JX n*N c a* a , 2 f *y /Q» V £T(^ ^ L (l) 

dxjdk T 4k T Q2 Uff[XJ,T) j2iir\k T J (l- 7 )r(2- 7 ) 

where we have defined the following Mellin-transforms 

#.,i(7) - / d 2 r(r 2 Q 2 y-~< J dz\^ T L (r, z; Q)\ 2 (10) 

which are given by 

/ 4( 7 ) \ 2a em N c v 2 1 r 2 (l+ 7 )r 2 (l-7)r 2 (2- 7 ) /(l+ 7 )(2- 7 )\ 

\<f>I(l) J " it ^ ? 4^7 r(2-2 7 )f(2 + 2 7 )(3-2 7 )i 2 7 (1 - 7) )' 1 ' 

Inserting formula J^Jl into 0) gives the forward-jet cross-section in the BFKL energy regime. One can easily show 
that the result is identical to the one obtained using kx— factorization pH l23|. The only undetermined parameter is 
a = a s N c /iT (with a s the strong coupling constant kept fixed) which appears in the exponential in formula ({§]). 



FtW) (9) 



C. The saturation regime 



Contrary to the BFKL case, the onium-onium cross-section in the saturation regime has not yet been computed 
from QCD. Studies are being carried out to identify the dominant terms in the multiple gluon exchanges [Til Izil 



25, 26] but the cross-section af at w ■. remains unknown. To take into account saturation effects, we are led to use a 

(??)(99) ... 1-3, . . 

phcnomcnological parametrization. We consider the following model introduced in |17| which is inspired by the GBW 
approach: 

=4^o(l-exp (-^^)) ■ (12) 
The dipole-dipole effective radius r 2 s (r,f) is defined through the two-gluon exchange: 

4™ 2 r 2 ff (r,f) ee of^, (r.r.O) = 4™ 2 min (r 2 f 2 ) {^S^^y} ■ ( 13 ) 

For the saturation radius we use the parametrization Rq(Y) = e~~s ^ Y ~ Y °> /Qo with Qq = 1 GeV. 
Let us express the cross-section in terms of a double Mellin-transform: 

< ) u^,Y ) =^ j ±lJ ^{-^m) 7 (4w) T - 9(7 ' T) (14) 

with 

gh,r)= [ du 2 [ du 2 u 2 "<- i u- 2T - 2 (l - e- r ifA u ' u A = 2F ( 7 ~ T ) {^(1,3+t-7,2t) + ^(1,3+t-7,2-2 7 )} 
Jo Jo v / l + r-7 

0<Re(r), Re{i), Re{<y-r) <1 (15) 

where the confluent hypergeometric function of Tricomi $(1, a, 6) can be expressed p7j in terms of incomplete Gamma 
functions. Inserting Ijl4|l in and (J2J), one obtains 



do- 1 T*r JX n 2 N c a 2 s <r 2 / d 7 2 2 . 7 . /" dr 2 2 4 T r 2 r(r) 

^yjfeff(xj,k T ) J — (4Q i? (y)) 7 ^ L (7) y — (4fc T i? (r)) -^j- 



ff(7»r) ■ (16) 



Inserting formula i|ltj|) into gives our parametrization of the forward-jet cross-section in the saturation regime. 
The parameters are A, Yq and the normalization <7q. 



D. Fixing the parameters 

The first sets of data published by the HI |14j and ZEUS collaborations regarded the measurement of da/dx. In 
previous studies, we fitted the BFKL-LL |l6j| and saturation parametrization Q on those data with the cut x<10~ 2 . 
Despite corresponding different energy regimes, in both cases we obtained good descriptions with x 2 values of about 
1. The obtained values of the parameters and the % 2 of the fits are given in Table I. 



fit 


parameters 


l/Ro(Y = 0) 


X 2 (/d.o.f.) 


BFKL-LL 
strong sat. 
weak sat. 


4a log(2)= 0.430 
A = 0.402 and Y = -0.82 
A = 0.370 and Y = 8.23 


1.18 Gev 
0.22 Gev 


12 (/13) 
6.8 (/ll) 
8.3 (/ll) 



TABLE I: Results of the BFKL and saturation fits to the first HERA forward-jet data. The saturation fits shows two independent 
solutions showing either strong or weak saturation parameters (see text). 

In the BFKL-LL case, the only parameter is a and the value obtained was 4alog(2) = 0.430. For the saturation 
fit, the two relevant parameters are A and Yq and the fit showed two \ 2 minima for (A = 0.402, Yq — —0.82) and (A = 
0.370, Yo = 8.23). We shall refer to the first (resp. second) solution as a strong (resp. weak) saturation parametrization. 
Indeed, the first saturation minimum corresponds to strong saturation effects as, for typical values of Y, the saturation 
scale 1/Ro is about 5 Gev which is the value of a typical kr- The second saturation minima corresponds to small 
saturation effects and rather describes BFKL physics. 

Along with formulae © and (|16J) . the values of the parameters given in Table I completely determine the BFKL-LL 
predictions and two parametrizations for the saturation model. We are now going to compare these with the very 
recent data without any adjustment of the parameters. This will provide a strong test of those small— x effects. 
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FIG. 2: The forward-jet cross-section da/dx. The points are measurement by the HI (left plot) and ZEUS (right plot) collabo- 
rations. The lines are comparisons with BFKL-LL predictions (full lines) and the two saturation parametrizations (dotted and 
dashed lines). In both cases, there is good agreement with the data. For comparison, fixed-order QCD predictions at NLO are 
also displayed. 



E. Comparison to the 2005 HERA data 
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FIG. 3: The forward-jet cross-sections da/dQ 2 (left plot) and da/dkr (right plot). The points are measurement by the ZEUS 
collaboration. The lines are comparisons with BFKL-LL predictions (full lines) and the two saturation parametrizations (dotted 
and dashed lines). For the two observables, there is good agreement with the data. 



We want to compare the cross-section Q obtained from the BFKL-LL prediction 10 and the saturation parametriza- 
tion with the new data coming from measurements performed at HERA [T^. |20| . On one side, our theoretical 
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FIG. 4: The forward-jet cross-section da / dxdQ 2 dk^ . The points are measurement by the HI collaboration. The lines are 
comparisons with BFKL-LL predictions (full lines) and the two saturation parametrizations (dotted and dashed lines for strong 
and weak saturation respectively). In the regime where on expects small— a; effects to be important (r = k^/Q 2 ~ 1), there is a 
good description of the data. In the regime where r ^> 1, the small— a; parametrizations do not reproduce the data as expected 
from the hierarchy of the hard scales. 



results are for the cross-section Q which is differential with respect to all four kinematic variables x, Q 2 , xj, and 
kx- On the other side, HERA data concern observables which are less differential: da/dx, da/dQ 2 , da/dk\, and 
da / dxdQ 2 dk 2 -, . Therefore, on the top of the Mellin integrations in 10 and (|16|) . one has to carry out a number of 
integrations over the kinematic variables which have to be done taking into account the kinematic cuts applied by 
the different experiments. A detailed description of how we performed those integrations in given in Appendix A 
and the resulting cross-sections that can be compared to the data are given in Appendix B. The method allows for a 
direct comparison of the data with theoretical predictions but it does not allow to control the overall normalization. 
In the following studies, we therefore compare only spectra and will not refer to normalizations anymore. As already 
mentioned, one does not adjust any of the parameters of Table I. 

Let us start with the observable da/dx which has been measured by both the HI and ZEUS collaborations and 
which now features lower values of x than the first measurements. The comparison is displayed in Fig. 2, the three 
small— a; parametrizations describe very well the data. One cannot really distinguish between the three curves, except 
at small values of x where one starts to see a difference: the BFKL curve is above the weak-saturation curve which 
is itself above the strong-saturation curve. However the main conclusion is that the data seem to feature the BFKL 
growth, when going to small values of x. For comparison, the fixed-order QCD predictions at NLO computed in 



with the DISENT Monte-Carlo program 28] are reproduced in Fig. 2. At the lowest values of x, they do not reproduce 
the data as they are about a factor 1.5 to 2.5 below depending on the experiment and the error bars. Even adding 
a resolved-photon component to the NLO predictions 20] does not pull them within the uncertainties, contrary to 
what happened for the previous data [2j|. This is an interesting difference with the forward-pion case y(j for which 
it seems that no higher-order effect other than a NLO resolved-photon contribution is needed. 

In Fig. 3 are represented the two other single differential cross-sections that we shall briefly discuss: da/dQ 2 and 
da/dkx measured by the ZEUS collaboration. One can see again that the three small— x parametrizations agree well 
with the data, it is a strong result that one is able to describe the Q 2 and kx spectra without any adjustment of the 
parameters as they were only fitted to describe the x dependence. 

We shall finally compare our predictions with the triple differential cross-section do j dxdQ 2 dk\ measured by the HI 
collaboration. The interesting part of this measurement is that it has been carried out in 9 different bins oir = k\jQ 2 
from 0.1 < r < 1.8 to 9.5 < r < 80. This allows to test the limits of our parametrizations which are supposed to be 
valid only when r ~ 1 as they do not take into account any transverse momentum ordering of the gluons emitted in 
rapidity between the forward jet and the photon. The comparisons with the data are shown on Fig. 4 and one sees 
the expected trend. The bins which have r ~ 1 are well described by the small— x parametrizations while the others 
are not: for the latter, we overshoot the data as the BFKL rise towards small values of x is too steep. Interestingly 
enough, the trend is reversed for QCD predictions at NLO: they describe better the data which feature large values 
of r. These observations favor the need of the BFKL resummation to describe the r ~ 1 data. The large— r bins also 
exhibit a limitation of the saturation model as one can see that the strong saturation parametrization lies above the 
weak saturation parametrization. Such a behavior indicates that the model should not be used when k\ 3> Q 2 ■ 

Let us comment further on the two saturation parametrizations. While the BFKL formula is a QCD prediction 
as it is computed from Feynman diagrams [l|, the saturation parametrization (|16J) is a phenomenological model. The 
fact that it describes well the data does not call for the same conclusions as in the BFKL case. It only exhibits that, 
as it is the case for a number of observables, data are compatible with saturation effects even at energies which do 
not require them. In other words, the forward-jet measurement at the present energies cannot distinguish between 
saturation and BFKL effects; one would start seeing a significant difference at higher energies. 

It is the purpose of the next section to look for such differences, by studying another process similar to forward-jets, 
namely Mueller-Navelet jets, at LHC energies. 



III. TOWARDS THE LHC: MUELLER-NAVELET JETS 



Mueller-Navelet jet production in a proton-proton collision is represented in Fig. 5 with the different kinematic 
variables. We denote y/S the total energy of the collision, k\ and &2 the transverse momenta of the two forward jets 
and x\ and X2 their longitudinal fraction of momentum with respect to the protons as indicated on the figure. In the 
following, we compute the Mueller-Navelet jet cross-section in the high-energy limit, recall the BFKL predictions and 
formulate our saturation model. We then display predictions for observables which can be measured at the LHC. 



Formulation 



As in the original paper Q, we consider the cross-section differential with respect to x% and X2 and integrated 
over the transverse momenta of the jets with k\ > Q\ and k% > Qi- Qi and Q2 represent then experimental kr— cuts. 
Considering the high energy limit, the QCD cross-section for Mueller-Navelet jet production reads [T^. l3lj: 



d<jP^ JXJ ir 2 N< 



dx\dx2 64 



r-OO r-OG 

feff(x 1 ,Ql)f eff (x 2 ,Ql) dr df QiJ 1 (Qir)Q 2 Ji(Q2f) o (gg)(gg) (r,f,Y) (17) 
Jo Jo 



with Y — \og(xiX2S/QiQ2) the rapidity assumed to be very large. 0( gg )(gg) (r, r, Y) is the cross-section in the collision 
of two gg dipoles of size r and f with rapidity total Y. As before / e // is the effective parton distribution function JSJ. 

Let us comment formula (|17H . As before, each forward jet involves perturbative values of transverse momenta and 
moderate values for x\ and x 2 . This explains the collinear factorization of the two functions f e ff, here we have taken 
the factorization scales to be Q\ and Q 2 . The remaining hard interaction is between two gg dipoles: as we have 
seen in the previous section, each of them describes a gluon emission at high energies. Formula l|17|) expresses the 
Mueller-Navelet jet observable in terms of the cross-section <7( g9 )( gs ) which contains the high-energy QCD dynamics. 
This is the similarity with the forward-jet case: the problem is also analogous to the one of onium-onium scattering. 





FIG. 5: Mueller-Navelet jet production in a proton-proton collision. The kinematic variables of the problem are displayed. S 
is the total energy squared, k\ and &2 are the transverse momenta of the jets and xi and xi are their longitudinal momentum 
fraction with respect to the incident protons. Y is the rapidity interval between the hard probes. 



Let us first consider the BFKL energy regime, the gg-gg dipole-dipole cross-section reads 
a BFKL (r = Y)- 2lrN - a2 s r 2 f rf 7 (r/rf 7 f a s N c \ 

which combined with (jl7[) gives 

dfJ BFKL _ n 3 N 3 a 2 r ^ (Q 1 /Q 2 ) 2 ^ / a s N c 



dx\dx2 SCfQI 



One can easily show that the result is identical to the one obtains using kx— factorization Q. As in the forward-jet 
case, the only undetermined parameter is a which appears in the exponential in formula 1)19(1 . We shall consider in 
this study the same value that was used for forward jets, that is a = 0.16. 

For the saturation paramctrization, we use the following gg-gg dipole-dipole cross-section: 

which up to the normalization is the same as ^12(1. The effective radius r e ff is defined by formula (|13J) and the 
saturation radius by Rq(Y) =e~3 ^ Y ~ Y °^/Qo with Qo = l GeV. Inserting (j20J into l(T7|. one obtains |3lj 

da sat n 3 N 3 a 2 a ( . f°° du f 2Q 1 Q 2 uRq(Y) 



^f eff (x 1 ,Q 2 1 )f eff (x 2 ,Q 2 2 )^l-2R 2 (Y)Q 1 Q 2 ^ 1+ f g ( u ) h 



d Xl dx 2 16C F "jjn^ujijjv"*,-**, ^ -w,-*^^ 1 + log(u) -y l+log(«) 

ex P ~- m — r-r R o( Y ) + ex P — m — r^ R a( Y ) 



l + log(u) u ' 7 V l + log(u) 



(21) 



In the following we consider only the strong saturation parametrization to display what could be the maximal expected 
effects at the LHC. The parameters are A = 0.402 and Yq = —0.82. The normalization oq is a priori not determined 
but we have fixed it so that at large momenta and small Y, one obtains the BFKL result. 




FIG. 6: The Mueller-Navelet jet cross-section do j dy\dyi as a function of j/i for different values of Q\. The kinematics of the 
other jet are fixed at Q2 = 30 GeV and j/2 = —4.5. The full lines are BFKL-LL predictions and the dashed lines are the saturation 
parametrization. 



B. Phenomenology 

We are going to study the dependence of the cross-sections Ijl9|) and H21JI as a function of the different kinematic 
variables x\, X2, Qi, and Qi- We want to consider large rapidities Y which implies very forward jets and therefore 
large values of x\ and X2- The well-known problem is that the cross-section is then damped by the parton distribution 
functions which at large x become very small. This prevents one to see the BFKL enhancement of the hard part of 
the cross-section with rapidity. 

A way out of this problem is to consider the following observables Ti-g/§ '■ 



- "l^T^, S)/^^(Qi,Q 2 , S) , (22) 



in other words, cross-section ratios for same jet kinematics and two different values of the total energy squared (S 
and S). The advantage of such observables is that they are independent of the parton densities and allow to study 
more quantitatively the influence of small— x effects [la UJl . ^ or instance the BFKL-LL prediction is (via a saddle 
point approximation) : 



g X 4alog(2) 



K s/ s -(j) ■ (23) 

The experimental verification of this at the Tevatron |32j | was not conclusive. The data were found above the prediction 
however it has been argued that the measurement was biased by the use of upper kr— cuts, the choice of equal 
lower kr— cuts, and hadronization corrections. The ratios (|22|l also display in a clear way the saturation effects |17ll31| 
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FIG. 7: Ratio of the saturation and BFKL-LL Mueller-Navelet jet cross-sections da / ' dy\dyi as a function of j/i for different 
values of Q\. The kinematics of the other jet are fixed at Q2 = 70 GeV and y2~— 3.5. 



which lead to ratios that, as a function of Qi and Q2, go from the value 1|23|) to 1 as the momentum cuts decrease 
into the saturation region (see Fig. 4 in Ref. 01 an d Fig. 3,4 in Ref. [3l|p. 

There is however an important experimental limitation to carry out the measurement H22JI : it would require to run 
the LHC at two different center-of-mass energies. If it turns out not to be possible, then one should settle for the 
cross-section da j dx\dxi- We shall now exhibit some of its characteristics, fixing the LHC center-of-mass energy at 
\/~S = 14 TeV. Also the absolute normalization is fixed to reproduce the Tevatron data at 1-8 TeV published 

in |32|. These data feature somewhat large error bars which leads to a significant uncertainty on the normalization 
for the LHC predictions. 

Let us introduce the rapidities of the two jets: 

»— < M > 

We first considered the case where one of the two jets has fixed kinematics Q2 = 30 GeV and j/2 = — 4.5. We looked at 
the dependence of the cross-section da jdy\dyi — x\Xi da j dx\dxi as a function of the other jet kinematic variables. 
In Fig. 6, we plotted the results for the BFKL-LL prediction l|19|) and the saturation parametrization pijl where the 
different plots feature the y\ dependence for different values of Q±. As expected, the cross-sections decrease quickly 
as yi gets large which corresponds to x± getting closer to one. For each value of Q\, one cannot really see a difference 
between the behaviors of the BFKL and saturation curves as a function of y±. However the relative normalization 
between the two curves is quite sensitive to the value of Q\. This is better exhibited on Fig. 7 where one displays the 
ratio of the saturation and BFKL results of Fig. 6. The ratio goes down to about 0.3 for Q\ = 10 GeV which represents 
a significant difference between the BFKL and saturation predictions. Note that this difference does not appear to 
be that large on Fig. 6 where the cross-sections are plotted. 

The second case we considered is the symmetric case Q = Qi—Q 2 and y = yi — ~y 2 which allows to go to bigger 
values of Y. We looked at the dependence of the cross-section da jdyidy^ as a function of Q and y. In Fig. 8, we plotted 
the results for the BFKL prediction (|19|) and the saturation parametrization i|21|) where the different plots feature 
the Q dependence for different values of y. In this case, the cross-section falls even faster when y gets big as both X\ 
and X2 get close to 1. Again, because of that, one does not see on the plot the difference between the BFKL and the 
saturation curves, yet it is still quite big as shown on Fig. 9 where we have displayed the ratio of the saturation and 
BFKL results of Fig. 8. For j/ = 5.5 and Q decreasing down to 10 GeV, the ratio goes down to about 0.4. 

We did not include in this study the weak saturation parametrization, the corresponding curves would lie in between 




FIG. 8: The Mueller-Navelet jet cross-section dajdyxdy^ as a function of Q = Q\ = Qi for different values of y = yi = — j/2- The 
full lines are BFKL-LL predictions and the dashed lines are the saturation parametrization. 



the BFKL and strong saturation curves which are displayed, and even is closer to the BFKL curve. There is a number 
of other plots one could study showing other dependences of da/dyidy2 but they are not needed for drawing our 
conclusions: testing BFKL effects and saturation effects with the observable dajdyxdy^ at the LHC will be a major 
experimental challenge as one will have to measure cross-sections with a great precision. We insist that this is due to 
the fact that the parton distribution functions at large x really damp the cross-section. Obtaining a high accuracy 
is not unfeasible because of the high luminosity at the LHC but this will require a very good understanding of the 
systematics errors. However we would like to emphasize the fact that better tests of small— a; effects could be realized 
with the measurement of the ratio 7Z>g/§, see formula (122(1 . 



IV. CONCLUSIONS 



Let us summarize the main results of the paper. The first part of the work was devoted to the study of the forward- 
jet measurement. We started by computing the QCD cross-section for forward-jet production in the high-energy 
(small— a;) limit. We recalled the BFKL-LL predictions JJjJ and also formulated the phenomenological model l|16|l that 
takes into account saturation effects. We then compared the BFKL and saturation-model predictions to the recent 
data from HERA for a number of observables: da/dx, da/dQ 2 , da/dh^, and da j dxdQ 2 dk\ . 

We obtained a very good agreement with the BFKL predictions and saturation parametrizations also show com- 
patibility with the data. Along with the fact that QCD at NLO predictions do not reproduce the small— a; data, this 
observation leads us to the conclusion that the present forward-jet data display the BFKL enhancement when going 
to small values of x. 

However, to make a definitive statement, one would have to make comparisons with BFKL predictions at next-to- 
leading-logarithmic (NLL) accuracy. The latter are under investigations and will hopefully be available soon, fn the 
mean time, let us discuss the expected qualitative impact of these BFKL-NLL corrections. Because we are describing 
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FIG. 9: Ratio of the saturation and BFKL-LL Mueller-Navelet jet cross-sections da / 'dy\dy^ as a function of Q = Qi = Q2 for 
different values of y = y\ = —y2. 



a kinematic regime in which fc|,~Q 2 , one can infer that they should be small. For instance, the contribution coming 
from the running of the coupling between those two scales would be unimportant. By comparison, the BKFL- 
NNL corrections seem to be very large for the proton structure function measurement J2J], in which case one easily 
understands why: the evolution takes place in a large range, from the soft proton scale up to the hard scale Q 2 . The 
situation is much different for forward jets. 

In the second part of the paper, we investigated small— a; effects for Mueller-Navelet jets in the LHC energy range, 
using the parameters that successfully describe forward-jets at HERA. We compared the BFKL-LL predictions Ijl9|l 
with those of the saturation model H21|l and concluded that the measurement of the simple cross-section will require 
a great precision to test the different scenarios. We argued that a better option to look for small— x effects was to 
measure the ratio of cross-sections (1221 which implies running the LHC at two different energies. 

On longer time scales, the international linear collider would give the opportunity to measure the virtual photon- 
virtual photon total cross-section at very high energies. This would also offer great possibilities |35l |36( for testing 
the BFKL enhancement and the saturation regime of QCD. 

Acknowledgments 

The authors would like to thank Robi Peschanski for useful comments and fruitful discussions. 

Appendix A: on the integration method 

To compare the forward-jet cross-section Q obtained from the BFKL-LL prediction and the saturation 
parametrization (|16|l with the data for observables which are less differential (da/dx, da/dQ 2 , da/dk 2 ^, and 
da 1 1 dxdQ 2 dk 2 ^) , one has to carry out a number of integrations over the kinematic variables. They have to be done 
while properly taking into account the kinematic cuts applied by the different experiments. This Appendix deals with 
these issues. 

Let us start from the quadruple differential cross-section da/dxdQ 2 dxjdk^, see formula JJJ. First one performs the 
Mellin integrations of @ and Then we choose the appropriate variables for the remaining integrations: to avoid 
numerical problems in the integral calculations, we chose variables which lead to the weakest possible dependence 
of the differential cross-section. We noticed that the best choice is 1/Q 2 , l/fej, log(l/xj), and \og(l/x). Since the 



experimental measurements are not differential with respect to xj, we carry out the integration 



d^a f xjS$a „ (l 

a log 



dxdQ 2 dk 2 ~, J dxdQ 2 dx j dk 2 -, \Xj 

With the BFKL-LL formula 10 , the convergence is fast enough so that one can perform all the remaining integrations 
to obtain any of the four observables mentioned above. With the saturation formula l|16(l . because of the extra Mellin 
integration and the time it takes to compute the functions 'J, performing all the remaining integrations would require 
important numerical work. We chose to use another method to obtain the cross-sections in that case. We shall 
describe it now with the example of da / dx. 

For a given value of x, the first step is to compute the differential cross-section 



dx J dxdQ 2 dk 2 ^ \Q 2 J v^t, 

for the BFKL case. The second step is to compute the bin center (kj. c , Qq) defined as follows: 

d {3) a BFKL d(J BFKL /dx 



dxdQ 2 dk 2 T {Q2c,k2rc) - JdQ 2 dk, 



The bin center is thus the point in the (fcy, Q 2 ) phase space where the differential cross-section in k 2 ^ and Q 2 is equal 
to the integral over the bin divided by the bin size (we will specify the integration limits later on). The third step is 
to obtain the cross section for the saturation case. We compute the cross-section at the bin center (fc 2 ^, Qfj)- 



da sat d (3) a s 



dx dxdQ 2 dk 2 

This procedure is valid if the bin center does not change much between the BFKL and saturated cross-sections. In 
other words, it means that the difference between the BFKL and saturated cross-sections is small. We saw in Section 
II-E that this is indeed the case. The method is easily adapted to the case of da/dQ 2 for which one finds a bin center 
(xc, fc| c ) for each value of Q 2 and to the case of da /dk 2 -, with a bin center (xc, Qc) f° r each value of kx- 

For the triple differential cross-section da / dxdQ 2 dk 2 ^ which is measured as a function of x, integrating over xj is 
not enough since one does not know the (Q 2 , fef.) bin-center. Instead, for a given value of x, one integrates also over 
Q 2 and k 2 -. and then divide the result by the Q 2 and k\~, bin sizes to obtain da / dxdQ 2 dk 2 -, . This is done for the BFKL 
case and one uses again the method described above to compute the cross-section in the saturation case. 

The other difficulty arises when setting the integration limits as one has to take into account the correlations 
between the kinematic variables (for instance y — Q 2 /sx < 1) and the cuts applied by the experiments (for instance 
cuts on the forward jet phase space). There are two ways to take these into account: either appropriately set the 
limits of integration while computing the integrals or evaluate later the phase space correction due to the experimental 
cuts. We are going to use both, since it is not possible to include all experimental cuts while computing the integrals. 

Table II is a list of the different set of cuts used by the HI and ZEUS experiments to carry out their measurements. 
For the ZEUS cuts, we only consider what they call the "forward-BFKL phase space" which corresponds the most 
to the Regge limit kinematics. Refering to this table, let us enumerate the integration limits which are used for the 
different integral calculations: 

• da/dx for HI: We integrate over 1/k 2 -, with the limit on k 2 -, defined by 0.5 < k^/Q 2 < 5 (this is an extra cut 
that HI applies to this measurement only), over l/Q 2 with Q 2 <sx, and over log(l/xj) with 0.035 < x ,/< 1 

• da/dx for ZEUS: We integrate over 1/k 2 -, with the limit on k 2 -, defined by 0.5 < k\jQ 2 < 2, over l/Q 2 with 
Q 2 <sx, and over log(l/xj) with 2<log(l/xj) < 3 

• da/dQ 2 for ZEUS: We integrate over 1/k 2 -, with the limit on k 2 -, defined by 0.5 < k\/Q 2 <2, over log(l/x) with 
x>Q 2 /s, and over log(l/xj) between 2 and 3 

• da/dk 2 -, for ZEUS: We integrate over log(l/x) with x > Q 2 /s, over l/Q 2 with the limits on Q 2 defined by 
0.5 < k\,/Q 2 < 2, and over log(l/xj) between 2 and 3 



da/dxdQ dk? for HI: The and l/Q limits of the integrals are defined by the bin values measured by the 
HI collaboration with also the kinematic constraint 0.1 <y = Q 2 /sx < 0.7. The log(l/xj) limits are obtained 
taking into account the cuts on the forward-jet angle which leads to 1.7354 <log(l/xj) < 2.7942. 



HI 


ZEUS 


E e > 10 GeV 
0.1 < y < 0.7 
10" 4 < x < 4.10" 3 
5 < Q 2 < 85 GeV 2 

k T > 3.5 GeV 
7 <9j < 20 degrees 
xj > 0.035 


E e > 10 GeV 
0.04 < y < 1. 

Q 2 > 25 GeV 2 
k T > 6 GeV 
2<rjj <3 

0.5 < fc^/Q 2 < 2 



TABLE II: ZEUS ( "forward-BFKL phase space") [Tflj ] and HI [2(ij ] cuts to define the forward-jet phase space. _E e is the energy 
of the outgoing electron and r\j = log(l/a;./) = — logtan(#j/2). The other kinematic variables have been defined in the text. 

The effects of the cuts defined in Table II which are not used above need to be computed using a toy Monte Carlo. 
They are modeled by bin-per-bin correction factors that multiply the cross-sections obtained as described above. 
This is how one proceeds: we generate flat distributions in the variables 1/fcfi, 1/Q 2 , log(l/a;j), and log(l/:r) using 
reference intervals which include the whole experimental phase-space (the azimuthal angle of the jet is not used in the 
generation since all the cross-section measurements are independent of that angle). In practice, we get the correction 
factors by counting the numbers of events which fullfil the experimental cuts for each a;— bin when we compute da/dx, 
each Q 2 — bin when we compute da/dQ 2 and so on. The correction factors are obtained by the ratio of the number of 
events which pass the experimental cuts and the kinematic constraints to the number of events which fullfil only the 
kinematic constraints, i.e. the so-called reference bin. Of course the experimental or kinematic cuts which have been 
applied already while computing the integrals are not applied in this study to avoid double counting effects. 

This method allows for a direct comparison of the data with theoretical predictions but it does not allow to control 
the overall normalization. This would require a full Monte-Carlo. Note that we did not use one in order to avoid any 
strong model dependence of the correction factors as they are only due to kinematic-cut effects. The derivation of 
the correction factors is independent of the theoretical input. They are given in Appendix B and they can be used to 
test any model suitable for the forward-jet cross-section, providing the same integration method as described above 
is used. 



Appendix B: tables with correction factors and resulting cross-sections 

In this section, we list the corrections factors that we obtained for the observables da/dx (HI and ZEUS), da/dQ 2 
(ZEUS), da/dk\ (ZEUS), and da / dxdQ 2 dk\ (HI). We also give the resulting cross-sections for the different points 
that we used to draw the curves on Fig. 2, Fig. 3, and Fig. 4. 



X 


factor 


bfkl-11 


weak sat. 


strong sat. 












0.00015 


0.24 


1200. 


1046. 


897. 












0.0005 


0.81 


805. 


785. 


722. 




X 


factor 


bfkl-11 


weak sat. 


strong sat. 


0.0010 


0.86 


371. 


365. 


360. 




0.00075 


0.13 


39.3 


34.5 


31.1 


0.0015 


0.80 


205. 


202. 


203. 




0.0017 


0.43 


62.5 


58.1 


56.0 


0.0020 


0.66 


117. 


114. 


116. 




0.004 


0.48 


25.4 


24.7 


24.6 


0.0025 


0.54 


72.0 


70.4 


71.7 




0.01 


0.34 


3.33 


3.38 


3.34 


0.0030 


0.45 


47.1 


45.9 


46.9 




0.025 


0.13 


0.106 


0.109 


0.106 


0.0035 


0.38 


32.3 


31.4 


32.1 












0.0040 


0.31 


22.5 


21.8 


22.3 













TABLE III: Correction factors due to experimental cuts and the resulting corrected cross-sections for da/dx in nb for BFKL-LL, 
weak saturation and strong saturation (see Fig. 2). Left Table: for HI cuts, right Table: for ZEUS cuts. 



Q 2 (GeV 2 ) 


factor 


bfkl-11 


weak sat. 


strong sat. 




k T (GeV) 


factor 


bfkl-11 


weak sat. 


strong sat. 


35. 


0.60 


6.75 


5.36 


3.81 




7. 


0.32 


71.9 


72.4 


71.6 


65. 


0.70 


1.40 


1.03 


0.843 




9. 


0.31 


21.1 


20.2 


20.4 


150. 


0.88 


0.184 


0.139 


0.137 




12. 


0.28 


6.23 


5.59 


5.63 


350. 


0.88 


0.0127 


0.0112 


0.0112 




17.5 


0.22 


1.01 


0.828 


0.828 


950. 


0.96 


0.000611 


0.000549 


0.000549 




25. 


0.15 


0.119 


0.0875 


0.0875 



TABLE IV: Correction factors due to experimental cuts and the resulting corrected cross-sections for BFKL-LL, weak saturation 
and strong saturation. Left Table: for da/dQ 2 in pb/GeV 2 , right Table: for da/dkr in pb/GeV (see Fig. 3). 



k\ (GeV 2 ) 


Q 2 (GeV 2 ) 


X 


factor 


bfkl-11 


weak sat. 


strong sat. 


12.25 < k% < 35 


5 < Q 2 < 10 


0.0002 


0.065 


8.28 


7.85 


4.40 






0.0004 


0.065 


3.86 


3.69 


2.53 






0.0006 


0.065 


2.04 


1.96 


1.51 






0.0008 


0.065 


0.705 


0.674 


0.591 






0.001 


0.065 


0.160 


0.151 


0.151 


12.25 < k% < 35 


10 < Q 2 < 20 


0.0002 


0.042 


0.571 


0.551 


0.332 






0.0004 


0.065 


1.15 


1.12 


0.826 






0.0006 


0.065 


0.752 


0.735 


0.612 






0.0008 


0.065 


0.544 


0.530 


0.473 






0.001 


0.065 


0.416 


0.406 


0.376 






0.0012 


0.065 


0.289 


0.280 


0.275 






0.0014 


0.065 


0.171 


0.166 


0.169 






0.0016 


0.065 


9.90c-2 


9.53e-2 


9.95e-2 






0.0018 


0.065 


5.70e-2 


5.07c-2 


5.34e-2 






0.002 


0.065 


2.20e-2 


2.11e-2 


2.25e-2 


12.25 < k% < 35 


20 < Q 2 < 85 


0.001 


0.060 


4.79e-2 


4.72e-2 


3.97e-2 






0.0015 


0.060 


3.50e-2 


3.45e-2 


3.09e-2 






0.002 


0.063 


2.68e-2 


2.63e-2 


2.46e-2 






0.0025 


0.065 


1.86e-2 


1.82e-2 


1.76e-2 






0.003 


0.065 


1.20e-2 


1.17c-2 


1.15e-2 






0.0035 


0.065 


8.08e-3 


7.90e-3 


7.90e-3 






0.004 


0.065 


5.60c-3 


5.47e-3 


5.52e-3 



TABLE V: Correction factors due to experimental cuts and the resulting corrected cross-sections for da/dxdk^dQ 2 in nb/GeV 4 
(bins with 12.25 <k?< 35 GeV 2 ) for BFKL-LL, weak saturation and strong saturation (see Fig.4). 
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TABLE VI: Correction factors due to experimental cuts and the resulting corrected cross-sections for da / ' dxdk^dQ 2 in nb/GeV 4 
(bins with 35<fcy <95 GeV 2 ) for BFKL-LL, weak saturation and strong saturation (see Fig.4). 
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TABLE VII: Correction factors due to experimental cuts and the resulting corrected cross-sections for da / dxdk^dQ 2 in nb/GeV 4 
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